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Abstract 

Based on a sample of 5.8 x 10^ J/i/j events taken with the BESII detector, the 
branching fractions of J/il^ — 2(7r'*"7r~)77 and J/'^ — 3(7r+7r~)77 are measured for 
the first time to be (2.26 ± 0.08 ± 0.27) x lO^^ and (7.24 ± 0.96 ± 1.11) x 10"^, 
respectively. 

PACS: 13.25.Gv, 14.40.Gx, 13.40.Hq 
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1 Introduction 



More than one hundred exclusive decay modes of the J/%1) have been re- 
ported since its discovery at Brookhaven [1] and SLAC [2] in 1974. According 
to Rcf. [3], direct hadronic, electromagnetic and radiative decays represent 
roughly 65%, 14%, and 7% of the total J/il) decay rate, respectively Up to 
now, only about half of all hadronic decays, 34.8%, have been measured in 
exclusive reactions. The sample of 58 million J/V' events, taken at BESII, pro- 
vides a chance to measure some of the missing hadronic decays. In this analysis, 
we report the first measurements of J/%1) — > 2(7r"'"7r~)r^ and J/tl) 3(7r"'"7r~)ry. 

The upgraded Beijing Spectrometer (BESII) detector located at the Beijing 
Electron-Positron Collider (BEPC) is a large sohd-angle magnetic spectrome- 
ter which is described in detail in Ref. [4]. The momentum of the charged par- 
ticle is determined by a 40-layer cylindrical main drift chamber (MDC) which 
has a momentum resolution of (Tp/p=1.78%\/l + {p in GeV/c). Particle 
identification is accomplished by specific ionization (dE/dx) measurements in 
the drift chamber and time-of-fiight (TOP) information in a barrel-like ar- 
ray of 48 scintillation counters. The dE/dx resolution is adE/dx = 8.0%; the 
TOP resolution for Bhabha events is cttof = 180 ps. Radially outside of the 
time-of-fiight counters is a 12-radiation-length barrel shower counter (BSC) 
comprised of gas tubes interleaved with lead sheets. The BSC measures the 
energy and direction of photons with resolutions of ue/E ~ 21%\/E (E in 
GeV), a^p = 7.9 mrad, and = 2.3 cm. The iron flux return of the magnet 
is instrumentd with three double layers of counters that are used to identify 
muons. 

A GEANT3 based Monte Carlo package (SIMBES) with detailed consideration 
of the detector performance is used. The consistency between data and Monte 
Carlo has been carefully checked in many high purity physics channels, and 
the agreement is reasonable. 



2 Analysis of J/V' — ^ 2(7r+7r )ri 

This decay is observed in the topology 7r"'"7r~7r~7r~77. Events with four charged 
tracks and at least two isolated photons are selected. The selection criteria for 
charged tracks and photons are described in detail in Ref. [5]. Each charged 
track must be well fitted to a helix, originating from the interaction region of 
Rxy < 0.02 m and \z\ < 0.2m, and have a polar angle, 9, in the range | cos^| 
< 0.8. 

Isolated photons are those that have energy deposited in the BSC greater 
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than 60 MeV, the angle between the direction at the first layer of the BSC 
and the developing direction of the cluster less than 30°, and the angle between 
photons and any charged tracks larger than 10°. To eliminate tracks from 7 
conversions, the minimum angle between any two oppositely-charged tracks is 
required to be greater than 10°. 

A 4C kinematic fit is performed under the hypothesis 7r^7r~7r^7r~77, and 
the chi-squared, x^^^-^+^-^^, is required to be less than 15. x^+^-^+^-^^ 
is also required to be less than the chi-squares for the K~^K~7r'^7r~^^ and 
7r+7r~7r'^7r~777 hypotheses. 




Fig. 1. The distribution of m^^ for candidate J/ip ^ r/2(7r+7r~) events. Dots with 
error bars are data, the histogram is the shape of the rj from Monte Carlo simulation, 
and the curve is background. 



After the above selection, the m^^ distribution is shown in Fig. 1, where a 
clear r] signal is observed. Monte Carlo simulation is used to estimate the 
background, and backgrounds from simulated channels are hsted in Table 1, 
where N^f is the number of events after event selection and N'^'^"^ is the back- 
ground normalized to 58 million J/ ip events. The sum of background events is 
18 events, which can be ignored. Another possible background channel is from 
J/V' — > 7?72(7r"'"7r~). No obvious rj signal is seen in the m^^ distribution from 
J/ip ^ ■yr)2{7T~^Tr~), as shown in Fig. 2(a). Therefore the background from this 
channel can also be ignored. 

Background from events with a Kg in the final state are estimated from rj 
sidebands. Fig. 2(b) shows the mass distribution of all 7i~^7i~ pairs for events 
with m-y-y in the rj region {\m^^ — 0.55| < 0.05 GeV/c^), and the shaded 
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Fig. 2. (a) Distribution of m-y-y from Monte Carlo simulation of J/tp "f2{TT~^TT~)ri. 
(b) Distribution of rriT^+T^- from the r] signal region (blank histogram) and r] side- 
bands (shaded histogram) for J — > 772(71+71") events. 



Table 1 

Background estimates for the J /ip r/2(7r+7r^) case. Nf^f are the number of events 
passing selection criteria; A^"'"'"^ is the background normalized to 58 million J /tp 
events. 



Channel 


MC sample 




jynorm 


J/ijj ^ 4>ri 


40,000 








J/ip prj 


40,000 








J/ip ^ 7r/7r"^7r~ 


70,000 


2 


3 


J/tp ^ ior]{i] 77) 


40,000 


11 


8 


J/^p — > ijjr( {r( tt^tt^t]) 


100,000 


135 


2 


J/ip ^ 4>ri'{r]' — > 7r^TT~r]) 


40,000 


2 







50,000 


6 


5 



histogram is for r] sidebands (0.45 GeV/c^ < m.y.y < 0.50 GeV/c^ and 0.60 
GeV/c^ < m^^ <0.65 GeV/c^). From Fig. 2, we conclude that the Kg signals 
are consistent with coming from background, as estimated from 7] sideband 
events. 



The m^^ distribution, shown in Fig. 1, is fitted with a Monte Carlo determined 
shape for the rj and a second order polynomial and yields 4839±158 J/ip ^ 
2{n^7i^)ri, Tj ^ 'j'j events. 
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3 Analysis of J/i/? — > 3(7r+7r )rj 



Events with six good charged tracks and at least two isolated photons are 
selected. The angle between two oppositely charged tracks is required to be 
greater than 10° to remove 7 conversions. A four constraint kinematic fit is 
made to the hypothesis J/i> ^ 3(7r+7r~)77, and xl+n-n+n-n+n-'yj required 
to be less than 15. 

After the above selection, the two photon invariant mass distribution is shown 
in Fig. 3; an 77 signal is evident. 




0.40 0.50 „0.60 0.70 

m(Yy) (GeV/c^) 

Fig. 3. Distribution of m^j for J/^ — 3(7r"''7r~)77 candidate events. Dots with error 
bars are data, the histogram is the rj shape determined by Monte Carlo simulation, 
and the curve is background. 

Monte Carlo simulation indicates that background from the decay modes 
listed in Table 1 can be ignored. Other possible backgrounds are from J/V' — 
^2{7r'^7r~)rj and J/t/j — > ^S{7r'^7r~)rj events. As described in Section 2, the m-y-y 
distribution from J/i/j ^ 72(7r"'"7r~)?7 shows no clear 77 peak. For J/i/j — >■ 
73(7r"'"7r^)?7, the m^^ distribution from Monte Carlo simulation, shown in 
Fig. 4(a), also does not show a peak in the rj region, so its contribution can 
also be ignored. 

The background with Kg final states can also be estimated, as was done 
previously, using the rj sidebands. Fig. 4(b) shows the tt+tt^ mass distribution. 
The full histogram is the 772,^+7^ distribution for events in the 77 signal region 
(|777^^ — 0.55| <0.05 GeV/c^), and the shaded histogram is for events from the 
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7] sidebands (0.45 GeV/c^ < m^^ < 0.50 GeV/c^ and 0.60 GeV/c^ < m^^ < 
0.65 GeV/c^). As in Section 2, we can conclude from Fig. 4(b) that most events 
with K'^ are associated with backgrounds which are measured by rj sidebands. 




m(Yy) (GeV/c^) m(ji^jt) (GeV/c^) 



Fig. 4. (a) The m-y-y distribution from Monte Carlo simulation of J ftp 73(7r+7r^)7/; 
(b) The distribution of 171^^+^- from ij signal events (histogram) and ij sidebands 
(shaded histogram) for J/ip ^ 773(7r+7r~) events. 

Fitting the 77 mass distribution in Fig. 3 with the Monte Carlo shape for 
1] and a second order polynomial, as described in Section 2, yields 616± 82 
events. 



4 Detection efRciency 

Initially events were generated according to uniform phase space. However, the 
cos 6' distribution of charged tracks in J/ip 2(7r+7r~)?7 was inconsistent with 
that from Monte Carlo simulation. Much better agreement is obtained when 
the angular distribution is generated according to 1+a cos^ 9, where a = 0.65± 
0.03 is obtained by fitting the cos6' distribution of charged tracks. Fig. 5(a) 
and Fig. 5(b) show the comparison of the cos 6 distributions for charged tracks 
and 1], respectively, with Monte Carlo simulated events with the charged tracks 
generated according to this distribution. Including the contribution from t] 
sidebands, the angular distributions are consistent. 

For J/ip 3(7r+7r^)?7, the detection efficiency is obtained from phase space 
events since the cos^ distributions of charged tracks and 77 are consistent 
with those from Monte Carlo simulation, as shown in Fig. 6(a) and Fig. 6(b), 
respectively. For the above two decay modes, the detection efficiencies are 
listed in Table 3. 
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Fig. 5. The cos 9 distribution of (a) 
charged tracks and (b) 77 particles, 
where the crosses arc data, the shaded 
histogram is from Monte Carlo simula- 
tion of J/ijj — 2{Tr'^Tr~)r), and the full 
histogram is the sum of sideband back- 
ground and Monte Carlo simulation. 



Fig. 6. The cos^ distribution of (a) 
charged tracks and (b) rj particles, 
where the crosses are data, the shaded 
histogram is from Monte Carlo simula- 
tion of J/ijj — 3(7r"'"7r~)j7, and the full 
histogram is the sum of sideband back- 
ground and Monte Carlo simulation. 



5 Systematic errors 



The systematic errors mainly come from the following sources: 

(1) The MDC tracking efficiency has been measured using channels like 
J/ip ^ AA and ip{2S) — >• 7:~^T:~J/ip, J/tl) ^ A*^A*~- It is found that the Monte 
Carlo simulation agrees with data within 1-2% for each charged track. There- 
fore, 8% and 12% are taken as the systematic errors in the tracking efficiencies 
for the 4-prong and 6-prong final states analyzed here. 

(2) The photon detection efficiency has been studied with different methods 
using J/ip ^ p^tt'^ events [7]; the difference between data and Monte Carlo 
simulation is about 2% for each photon. In this analysis, 4% is taken as the 
systematic error for the 77 decaying into two photons. 

(3) The kinematic fit is useful to reduce background. The systematic error 
from the kinematic fit is studied using the clean channel J/ip ^ 7r+7r~7r'^, as 
described in Ref. [5]. The efficiency difference of the kinematic fit between 
data and Monte Carlo simulation is about 4%. With the same method, the 
decay modes J/ip ^ 2(7r+7r^)7r° and J/ip 3(7r+7r^)7r° are also analyzed. 
The efficiency difference of kinematic fit between data and Monte Carlo are 
4.3% and 5.5% respectively Since J/i/j 2{tt+tt-)tt^ and J/i/j 3(7r+7r-)7r° 
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are similar to the two channels analyzed in this paper, 4.3% and 5.5% are 
taken as the systematic error of the kinematic fit. 

(4) Other possible decay modes which may contribute to r] signals 
have been studied, and the background from them can be ignored. The error 
from the background under the 7] peak is included in the fitting error. The 
uncertainties of the background shape in the two channels are estimated to 
be about 3.4% by changing the order of the polynomial. Possible background 
from the continuum events [8] is estimated using data at i/s = 3.07 GeV/c. 
After applying the same selection criteria as above, no significant rj signal is 
observed. Therefore, the background from this source is also negligible. From 
the above analysis, the background uncertainty for the two decay modes is 
less than 5%, which is taken as the background systematic error. 

(5) The branching fraction of 77 —> 77 is (39.43 ±0.26)% [6]. The error is also 
taken as a systematic error. 

(6) The number of J/'4> events is (57.70±2.72) xlO^, determined from inclusive 
4-prong hadrons [9]. The uncertainty, 4.7%, is also a systematic error. Table 
2 lists the systematic errors from all sources. 

Table 2 

Summary of systematic errors (%) 



Source 


2{tt^tt )ri 


3(7r"'"7r )r) 


MDC Tracking 


8 


12 


Photon efHeciency 


4 


4 


Kinematic fit 


4.3 


5.5 


Background 


5 


5 




0.7 


0.7 


Number oi J/ip evetns 


4.7 


4.7 


Total 


12.1 


15.4 



6 Results 

The branching fractions are calculated with the following relation: 

- n{n-'n-)r,) = ^"'^ , (1) 

s ■ B[ri 77) • Njj^ 

where n is 2 or 3, Nobs is the observed events, e is the detection efficiency, 
B{ri 77) is the branching fraction of r/ — * 77, and Nj/^ the total number 
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of J/V' events. 

Table 3 summarizes the quantities used in the calculation of the two branching 
fractions and the final results, including systematic errors. 

Table 3 

Numbers used and branching fractions measured. 



Decay Modes 


Nobs 


s{%) 


Branching Fraction 


J/V' 2{Tr+Tr~)ri 


4839±158 


9.43±0.10 


(2.26 ± 0.08 ± 0.27) x lO'^ 


J/V 3(7r+7r-)r? 


616±82 


3.74±0.06 


(7.24 ±0.96 ±1.11) X 10-^ 



7 Summary 

In this paper, the decays of J/V' ^ 2(7r"'"7r~)?7 and J/ip ^ 2(7r+7r~)?7 are stud- 
ied with the BESll 5.8 x 10^ J/tjj event sample and their branching fractions 
are measured for the first time to be: 

S( J/V' ^ 2{n+TT-)ri) = (2.26 ± 0.08 ± 0.27) x lO'^ 

S( J/V' ^ 3(7r+7r-)7/) = (7.24 ± 0.96 ± 1.11) x 10"^ 

Comparing with other branching fractions of J/ ip decaying into stable hadrons, 
the branching fractions of J/ip — >• 2(7r"'"7r~)r; and J/ip — >• 3(7r"'"7r~)?7 are not 
large. 
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